We intend to study how a blending of a tactical (i.e. direct response to resource availability or perception) and a strategic (i.e. memory-driven and forward-thinking) behavior can help a forager navigate a dynamic and heterogeneous resource environment, with a particular emphasis on the aquisition and adaptabilty of migratory behavior.
Out model is an advection diffusion with two advective terms corresponding to a direct response to the environment and to memory. Importantly, the environment is approximately periodic, i.e. one in which the distribution of the resource at a given time is approximately equal to that resource distribution at a fixed previous time period, i.e. \(h(t) \approx h(t - \tau)\). This models annual variations, and the period corresponds to one year. The tactical process is straightforward taxis up a resource gradient, while the collective memory corresponds to a tendency of the organisms to repeat the behavior of the population in previous years. Thus:
\[{\frac{\partial u(x,t)}{\partial t}} = -\varepsilon {\frac{\partial^2 u(x,t)}{\partial x^2}} + \alpha \frac{\partial h(x,t)}{\partial x} + \sum_{i = 1}^n \beta_i \frac{\partial u(x, t - i\tau)}{\partial x}\]
where \(u\) represents a population distributed in time and space, \(\varepsilon\) is a rate of randomness (diffusion), \(\alpha\) is the strength of advection in the direction of the resource \(h(x,t)\) gradient), and \(\beta_i\) is the strength of advection up the gradient of the population distribution in an \(i\)’th earlier seasonal cycle (year), up to \(n \in N\). The form of \(\beta_i\) should be simple, e.g., if the collective memory reaches only one year back then \(n = 1\) and \(\beta\) is a single constant. Alternatively, a form of memory might discounts older “traditions” via \(\beta_i = \beta \exp(-i/\gamma)\), where \(\gamma\) represents a generational scale of memory.
For reference and replicability, I include the incrementally constructed R code, which relies heavily on the deSolve (Soetaert et al. 2010) and ReacTran (Soetaert and Filip, 2020) packages, which efficiently implement transport ODE and PDE systems in R.
As a simple illustration of a model run, below the ultimate evolution of a forager that begins in the left part of space and uses a combination of the previous year’s location (to the right) and the current resource distribution, and ultimately “splits up” into being attracted to the resource in the middle and its previous year’s distribution, with a somewhat stronger taxis towards the resource than the memory:
This example is completely static and symmetric around the two drivers. Note that all the populations (previous, initial, final) are normalizes (integrate to 1) and that the space is defined between 0 and 100.
In the following example, we let the population be “initialized” via a sinusoidal migration pattern, while the resource is fixed in the middle of the space (as above). The memory is a one-year memory. In this context, migration is not, in fact, optimal. But the important piece is to see how the reponse mixes attraction to the resource with attraction to the previous year’s behavior. Note that the year here (for simplicity, speed and symmetry), is 100 days long.
The following examples illustrate how an initially migrating population evolves, even as it is attracted to a resource that is distributed in two different locations. In these examples, the first year (year 0) is a sinusoidal migration. We run the simulations for 3 subsequent years.
The population roughly echoes its initial track, but degrades from year to year.
In the following simulations, we add a resource distributed as follows in two concentrations:
Here, the population disregards its memory and only pursues a tactic stategy. Within a year, the population has gathered around the resource.
Finally, we combine both resource following and memory:
Here, the strong cultural memory is slowly blended with a strategy to simply concentrate in one of the two resource hot-spots.
There is obviously nothing particularly adaptive about migrating when resources are fixed. In the following scenaries, we explore a seasonal resource. We developed (with some fuss!) a seasonal resource function with the following properties:
To generate a resource with these properties, we distributing the resource in space as a beta distribution, where the two shape and scale parameters vary sinusoidally in such a way as to fulfill the criteria above. Thus:
\[R(x,t, \theta) = \chi B(x/\chi, a(t, \theta), b(t, \theta))\]
where \(\chi\) is the maximum value (domain) of \(x\), \(B(x, a, b)\) is the beta distribution, \(\theta\) represents the set of parameters \(t_r, x_r, \sigma_t, \sigma_x\), and the two shape parameters are given by:
\[a(t) = \frac{m}{s^2}( s^2 + m - m^2)\] \[b(t, x', \sigma') = (m-1)\left(1 + \frac{m}{s}(m-1)\right) \]
where \(m(t)\) and \(s(t)\) describe the dynamic mean and variance of the resource peak. These equations are solutions to the mean and variance of the beta distribution, \(\mu = \alpha/(\alpha + \beta)\), \(\sigma^2 = \frac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}\).
The means and variances themselves are Gaussian pulses, with the mean peaking at \(x_r\) at time \(t_r\) with standard deviation \(\sigma_t\) and at \(\chi - x_t\) at time \(\tau - t_r\) and the standard deviation pulsing from \(\chi/\sqrt{12}\) (corresponding to a uniform distribution) at times 0, \(\tau/2\) and \(\tau\) down to \(\sigma_x\) at \(t_r\) and \(\tau - t_r\), with standard deviation (in time) \(\sigma_t\).
All of these details boil down to the ability to generate various seasonal scenarios.
In these pulsed resources, we might assume that a combination of memory and tactic responses would yield more optimal behavion, and might observe some shift in the migratory strategy.
## running year 1
## running year 2
## running year 3
## running year 4
## running year 5
## running year 6
This model does a surprisingly good job aggregating at the appropriate resource, but appears to only locally aggregate, not to migrate per se.
## running year 1
## running year 2
## running year 3
## running year 4
## running year 5
## running year 6
It is not clear that the hybrid model recovers a true migration, but it seems to retain somewhat more concentration along the migratory corridors.
## running year 1
## running year 2
## running year 3
## running year 4
## running year 5
## running year 6
Here we compare a memory with and without memory exploiting the “narrow - long” resource combination, and the tendency to
Notable concentrations at the resource peaks, but a lot of diffuse behavior in the interim seasons
## running year 1
## running year 2
## running year 3
## running year 4
## running year 5
## running year 6
Memory appears to reinforce the aggregations at the (narrow) resource peaks
## running year 1
## running year 2
## running year 3
## running year 4
## running year 5
## running year 6
Karline Soetaert, Thomas Petzoldt, R. Woodrow Setzer (2010). Solving Differential Equations in R: Package deSolve. Journal of Statistical Software, 33(9), 1–25. URL http://www.jstatsoft.org/v33/i09/ DOI 10.18637/jss.v033.i09
Soetaert, Karline and Meysman, Filip, 2012. Reactive transport in aquatic ecosystems: Rapid model prototyping in the open source software R Environmental Modelling & Software, 32,